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RATES OF CONVERGENCE TO SCALING PROFILES IN A 
SUBMONOLAYER DEPOSITION MODEL AND THE 
PRESERVATION OF MEMORY OF THE INITIAL CONDITION 


FERNANDO P. DA COSTAt§, JOAO T. PINTOt§, AND RAFAEL SASPORTESt§ 

Abstract. We establish rates of convergence of solutions to scaling (or similarity) profiles in a 
coagulation type system modelling submonolayer deposition. We prove that, although all memory 
of the initial condition is lost in the similarity limit, information about the large cluster tail of the 
initial condition is preserved in the rate of approach to the similarity profile. The proof relies in a 
change of variables that allows for the decoupling of the original infinite system of ordinary differen¬ 
tial equations into a closed two-dimensional nonlinear system for the monomer-bulk dynamics and 
a lower triangular infinite dimensional linear one for the cluster dynamics. The detailed knowledge 
of the long time monomer concentration, which was obtained earlier by Costin et al. in [7] using 
asymptotic methods and is rederived here by center manifold arguments, is then used for the asymp¬ 
totic evaluation of an integral representation formula for the concentration of j'-clusters. The use 
of higher order expressions, both for the Stirling expansion and for the monomer evolution at large 
times allow us to obtain, not only the similarity limit, but also the rate at which it is approached. 

Key words. Dynamics of ODEs, coagulation processes, convergence to scaling behaviour, 
asymptotic evaluation of integrals, submonolayer deposition model 
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1. Introduction. Submonolayer deposition is the process of particle deposition 
onto a surface such that the deposited particles can diffuse and coagulate to form 
clusters, and those clusters are so sparse as not to cover the full original surface. 
The theoretical modelling of this process can be done using a variety of approaches 
[10]. Our work is a continuation of a number of recent mathematical studies of the 
dynamics of deposition using a mean-field approach [4, 5, 6, 7]. 

Denoting by Cj = Cj(t) the concentration number of clusters of size j, or j-clusters, 
at time t, on the deposition surface, assuming this surface to be bombarded at a 
constant rate a > 0 by 1-clusters (also called monomers)^ and considering the clusters’ 
concentrations so small that cluster-cluster reactions can be disregarded, the cluster 
dynamics on the surface can be modelled by a Smoluchowski’s coagulation system 
with Becker-Doring like coagulation kernel, meaning that the only allowed coagulation 
reactions are those in which a monomer takes part, namely (1) -|- (j) —> (j -I- 1), with 
rate coefficients oi.j > 0, and j € {1, 2,...}. 

An assumption that is relevant in some applications is the existence of a critical 
cluster size n below which clusters are not stable and do not occur in the system 
in any signihcant amount in the time scale of the coagulation reactions. There are 
several ways to model this assumption (see, e.g., [7] and references therein). In [7] it 
was modelled by considering that no clusters of size larger than 1 and smaller than 
or equal to n — 1 can exist, and thus the smaller cluster that is not a monomer has 
size n and is formed when n monomers come together and react into an n-cluster (if, 
as in Monte Carlo simulations, we consider the monomers sitting in the vertices of a 
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lattice, we could have these “multiple” collisions simply by having n — 1 monomers 
surrounding an empty site —as in the four nearest neighbours in a square lattice— 
which is suddenly bombarded by an aditional monomer, creating the n-cluster in a 
single n body reaction). 

If one considers that the coagulation rates are independent of the cluster sizes 
undergoing reaction (aij = 1, say), the mean-field model for submonolayers deposi¬ 
tion with a critical cluster size n becomes the following infinite system of ordinary 
differential equations: 


( 1 . 1 ) 


Cl = a — nci — Cj 

j^n 

Cn — C1C7.J, 

Cj = CiCj-i — ciCj, j > n + 1 . 


System (1.1) with n = 2 (i.e., with no unstable clusters) was considered in [ 4 , 5, 6], 
where the long time behaviour of solutions and the approach to a similarity profile 
was studied; in [ 4 , 6] it was even considered the case of time dependent monomer 
input a = a(t) of power law type. 

In [7] system (1.1) was considered with general n> 2 , and the results in [5] were 
correspondingly generalized to the following 

Theorem 1 ([7, Theorem 5]). Let (cj) be any solution to (LI). Consider the 
new time scale T{t) := ci{s)ds, and let Cj(T) := Cj(t(T)). Then, 


lim 

j, T —>■ -boo 
V = j /t fixed 
77 ^ 1 



(n —1) /n^ 

Cj{r) 




(1 — 77) , if 0 < 77 < 1 

0, if 77 > 1. 


This result means that, in the self-similar (or similarity) variable rj = j /r, the 
scaled solutions approach a universal profile $1(77) in the similarity limit j, r —>• 00 
with 77 constant. 

The purpose of the present paper is to study the rate at which the limit function 
$1(77) in Theorem 1 is approached. This is a kind of problem that is rather natural in 
long time dynamics, with the paradigmatic example being the linearization method 
which, when applicable, provides an exponential estimate for the rate of convergence 
of solutions to the limit (an equilibrium, say). 

In the context of Smoluchowski’s coagulation equations and analogous systems 
(such as coagulation-fragmentation, and Bercker-Doring systems) the proofs of ex¬ 
istence of, and convergence to, self-similar solutions are already so demanding that 
the study of the rate at which this convergence takes place is, so far, typically out of 
reach. 

To the best of our knowledge, so far only the two papers [ 2 , 11 ] considered the 
problem of the rate of converge to self-similar behaviour in Smoluchowski’s systems. 
In [2] the result was obtained for the constant kernel case (in continuous variables) 
ax,y = 1 using the approach of linearizing the coagulation equation in self-similar 
variables about the self-similar solution, and proving that the linearized operator has 
a spectral gap in an appropriate scale of weighted Sobolev spaces. In [ 11 ] the case of 
so called solvable kernels (i.e. = 1, x + y, xy) was considered using an approach 
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analogue to the Berry-Esseen theorem in classical probability theory. Both approaches 
require rather delicate and difficult analysis. 

In the present paper we are able to approach this rate of convergence problem 
using an essentially simpler approach due to the special nature of the coagulation 
reactions in (1.1). Exploring the approach which was already used to good effect in 
the proof of Theorem 1 in [5, 7], we can establish the main theorem of this paper: 

Theorem 2. Let (cj) he any non-negative solution of (1.1) with initial data 
satisfying 3pi,p2 > > 1 : Vj > n,j^Cj{0) C [pi,P 2 ]- Assume the notation 

introduced in Theorem 1. Then, as j, r —>■ +oo, with rj = j /t ^ 1 fixed, 


n — 1 



( 1 . 2 ) 



(1.3) 


(1.4) 


What we find interesting in this result is that, although Theorem 1 states that, in 
the similarity variable rj, solutions to ( 1 . 1 ) approach a universal profile independent 
of the initial condition (indeed this universal behaviour is what is physically relevant 
in this type of enquiry), the rate at which this limit behaviour is approached still 
preserves some information about the initial condition, namely the rate of decay of 
the initial datum for large cluster sizes, p, can still be computed from the observation 
of the decay rate of scaled solutions to $ 1 ( 77 ), for 77 > 1 . 

2. Preliminaries. Our approach to the study of system (1.1) follows the one 
used in [5, 7] and consists of the exploration of the following two observations: 

(i) First note that the equation for ci depends only on ci and on the “bulk” 

quantity y{t) := which (formally) satisfies the differential equation y = cf. 

Thus, the definition of this bulk variable allow us to decouple the resulting infinite 
dimensional system into a closed two-dimensional system for the monomer-bulk vari¬ 
ables (ci, y) , from which we get all the needed information about the behaviour of 

Cl. 

(ii) Secondly, the remaining equations for cj, with j > n, depend only on those 
same variables cj, and on ci. However, the way they depend on ci is such that, by an 
appropriate change of the time variable, the system is transformed into a linear lower 
triangular infinite system of ordinary diferential equations, which can be recursively 
solved, in terms of ci, using the variation of constants formula. 

Let us look more closely at the details of (ii) first: writing the second and third 
equations of ( 1 . 1 ) in the form 



it is natural to introduce a new time scale 


( 2 . 2 ) 



along with scaled variables 


(2.3) 


Ci(t) := Ci(t(r)) and Cj(r) := Cj(t(r)), j > n1 
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where t(r) is the inverse function of t(<). With this new time (2.1) reads 


(2.4) 


f •~ / ~n—1 •~ 

J — Cl Cji 

\ Cj' = Cj-i — Cj, j > n + 1, 


where (•)' = d/dr. The equation for c„ can readily be solved in terms of ci, 

Cn{r) = e~'^Cn{0) + f = e~'^Cn{0) + f - s)ds, 

Jo Jo 

and then we can use this solution to solve the c„+i equation, and afterwards solve 
recursively for j > n + 2. The following expression for Cj{T) is obtained in [7], 
generalizing an analogous result first proved in [5] for the special case n = 2: 


i{t-) = 


rJ-k 


3 

^ (j - k)\ 

k=n ' 


Cfe(O) 


(i - n) 


d 


(ci(r —s))” ^ "e ®ds, 


(2.5) 


Vj > n. 


Now, the information about ci needed to study (2.5) has to be extracted from 
the two-dimensional system for (ci, y) refered to in observation (i). Let us see this 
with a bit more detail. Denoting x{t) := ci(t) and y{t) := on.e can easily 

conclude that these quantities (formally) satisfy the two-dimensional system 

(2.6) { x = a -nx--xy 

y y = X . 


It can be proved [7, Theorem 1] that if X]^nCj (0) < oo, then any solution of 
system (2.4)-(2.6) is also a solution of system (1.1). 

Clearly, the monomer-bulk system (2.6) form a closed two-dimensional system 
the solutions of which can be studied independently of what happens to the infinite 
system (2.4). In particular, the needed information about Ci in (2.5) will be obtained 
from (2.6). 

3. Center manifold analysis of the monomer bulk system. In this section 
we will study the long time behaviour of solutions to (2.6). The following result, 
proved in [7], holds: 

Proposition 3. [7, Proposition 1] Let {x{-),y{-)) be any solution of (2.6) with 
non-negative initial eondition. Then: 

(i) (x(t), y(t)) is positive when t > 0. 

(ii) {x(t),y{t)) exists for all t > 0. 

(iii) x{t) —>-0, as t ^ -|-oo. 

(iv) y{t) —>■ -boo, as t ^ -boo. 

(v) Of — cry —>■ 0, as t ^ -boo. 

The next result establishes the rate of decay of x{t) as t —>■ oo. 

Theorem 4. [7, Corrected Eg. (25)] Let {x{-),y{-)) be any non-negative solution 
of (2.6). Then, 


1 n(n — 1) 1 , 

t-— + + o(r 


as t ^ -boo. 


(3.1) x{t) 


n -b 1 


n -b 1 
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Remark 1. In [7] the coefficient of the t~^ term is slightly different from that in 
(3.1), in what we believe to be a minor mistake. 

Proof, (of Theorem 4). The proof of this result in [7] was done using asymptotic 
methods. Here, we follow the idea of the proof of the case n = 2 in [5, Proposition 
3.2], which uses center manifold methods to obtain the first term in the right-hand 
side of (3.1). The main difference relative to [5] being that now we need to use a 
higher order approximation to the center manifold in order to obtain the expansion 
with an additional term, (3.1). However, somewhat surprisingly, this has to be done 
in a two step process, as explained further down. 

Let us define a new variable v{t) := a — x{f)y{f), so that {x{t),v{t)) —>■ (0,0) as 
t — >■ -Loo. 

By changing variables (x, y) i—>• (x, v) system (2.6) becomes 


(3.2) 


X = V — nx 
av 

V = - X 

X 


,n+l 


V 

X 


-L anx" ^ — Twx" ^. 


Suppose x(0) 0. Otherwise, by Proposition 3 (i) and (ii) we know that x(t) > 0 for 

all t > 0, and so just redefine time so that x(0) becomes positive. Now change the 
time scale 

(3.3) t^C = Cit)-=[ -^ds, 

and define (x(C), v{())) := (x(t(C)), v{t{C))) , where t{()) is the inverse function of ({t). 
By Proposition 3(iii), we have x{t) —0 as t — >■ -Loo, and so also ( -Loo ast ^ -Loo. 
With the new time scale system (3.2) becomes 

, . ( x' = vx — nx"'+^ 

( x' = —av — +iP + anx" — nvxP, 


where (•)' = d/d(, or in matrix form 


X 

/ 

■ 0 

0 


X 


vx — 

V 

— 

0 

—a 


V 

+ 

—anfP — nvx^ 


Since, by Proposition 3 (iii)-(v), we know that when t —> -Loo all non-negative 
solutions {x,y) of system (2.6) converge to (0,oo) and u —>■ 0, the transformation 
above imply that the corresponding solutions (x, u) of (3.5) converge to (0,0) when 
C —>■ +00. 

Using standard results in center manifold theory [3, chap. 2] it is straightforward 
to conclude that (3.5) has a center manifold in a neighbourhood of the origin, that 
locally exponentially attracts all orbits, and is the graph of a function v = (/)„(x), 
where 

(3.6) (/)„(J) = nx’^ - -X^+^ + _ ff+±~2n+4 ^ 0(j3n+2)^ 

a a‘‘ 

In order to get an expression for the behaviour of solutions x(C) for large values 
of C, from which (3.1) can be deduced, we need to proceed in two steps: first, using 
the first two terms in the expression of the center manifold 4>n and the knowledge that 
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x{() —>■ 0 as C —>■ +00, we obtain a first, “lower order”, expression for the long-time 
behaviour of x{(^). Then, using the next term in the right-hand side of (3.6) and this 
“lower order” information, we obtain a “higher order” correction to the long-time 
behaviour of x{Q which will be sufficient to prove (3.1). 

Remark 2. It is worth noting that the two stage process just described cannot 
be abbreviated to a single step, taking from the start the first three terms in fin and 
the knowledge that ir 0 as C —>• -l-oo. Doing this one would arrive at (3.1) with 
the equality substituted by an asymptotic equality, meaning that the long-time limit 
of the ratio of the left-hand side of (3.1) by its right-hand side is equal to 1. Note 
that this is a weaker result than the one stated in (3.1), and, in fact, it is not strong 
enough to prove our main result. Theorem 2. 

Let us start to implement the idea described above. Considering the approxima¬ 
tion to the center manifold given by taking only the fist two terms in the right-hand 
side of (3.6), 

finifi) = -p C>(^2’^+2), 

Q 

the dynamics on the center manifold is given by 

j; = -isy3^0(5?2n+3) as C^+OO, 
a 

which can be written as 

(3.7) + asC^+oo. 

Knowing that x{() —0 as C —-l-oo, the right-hand side of this equation converges to 
1 and we conclude that, for all e: > 0, there exists a T > (o such that, for all ( > T, 
the following inequalities hold 

d f at \ 


Integrating these differential inequalities between T and ( we get 


(3.8) 


(1-£)(C-T)< 


a 


{n + 2)xc (n + 2)x 




where xt = Xc{T). Dividing (3.8) by ( and taking liminf and limsup we obtain 

C-S-+00 C->-l-oo 


1 — £ < lim inf ■ 


a 


C-s-i-oo (^n -\-2)( X, 


•n+2 


< lim sup ■ 

C—>-|-oo 


a 


■2)(x 




^ 1 + iS, 


which, due to the arbitrariness of e, means that lim^^+oo C *c(C)”~'"^ = ^-^d 

thus 

(3.9) ^ C-^+oo. 

Let us now consider a better approximation to the center manifold, obtained by 
considering a further term in the right-hand side of (3.6), namely, 

fin{fi) = nJ" - - -k C>(j2"+4). 

a 
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The dynamics on the center manifold is now given by 

(3.10) x', = 7 + 0(£g"+^) as C^+oo. 


Writing this differential equation as 


—ax'c ' ” + n{n — l)a; 
and observing that 


:n +3 ; 7^ 1 \~ 2 n +3 “ 1 + C’(a;”+^), as C + 00 , 




ds = 


a 


+ n(n—l)s^"+3 (n + 2)s"+^ 


n(n — 1) 
2^2 


■^n(s), 


where 


V'ti(s) := n{n - ^) J 


c,n—3 


-ds. 


5 — 


i,(n—1) 


we conclude that (3.10) can be written as 

i {{n + 2)x-+^ + ^ ^ ^ ^ ^ +'^- 

As in the previous approximation (3.7), we now need to estimate the right-hand 
side of (3.11), but instead of using only the information that a;c(C) —>■ 0 as C —>■ -Loo, 
we shall use (3.9). 

Observing that there exist constants K* > AT* such that the right-hand side of 
(3.11) can be bounded by 1 -I- Ar*x"+^ < 1 + a;"^^0(l) < 1 + we can write, 

for all C > T, 


(3.12) i + 


(n -I- 2)x, 


n+2 


+ + Mxc )) < 1 + K*x:+^ 


Integrating these differential inequalities between T and C > T and denoting Xc(T) 
by xtj we get 


C - T + iC, J((+2(s)ds < 


(3.13) < ( 
and 


a _|_ n(n—l) 




a I n{n-l) 

(n+2)x5,+^ 


■1pn{xT)^ , 


C-T + K* x^+^{s)ds > 

- {in+2^S7+‘^ + +'0n(^c)) ” ( 


(n+2“£”+^ + + Mxt)) ■ 


We now use (3.9) to estimate the integral of : Let e > 0 be fixed arbitrarily, and, 
if necessary, redefine T such that C^”+^((() € [1 — e, 1 -I- e], for all C > T. Thus 


n + 2 


(1 — e)(logC — logT) < J a;"+^(s)cis < 


71 -|- 2 


(l-b£)(logC-log T). 
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On the other hand, to estimate 4’n{Xc{C)) as C —>■ +oo first observe that, defining 

y ■= s, we can write ipniv) = -a For n > 3 the 

explicit expression of this last integral is known from [ 8 , Integrals 2.146-3, 2.146-4] 
and it is easily seen to be bounded when y —>■ 0. The case n = 2 is a little trickier: we 
can easily compute 

r s-1 4 4 

tjj 2 {s )=2 / -5 -^ds =-log =-logs+ 0 ( 1 ), as s 0 , 

d s - 2 « - 2^2 a 

y OL 

and to evaluate tpnixdC)) when ( +00 we now need to use (3.9): in the case 

n = 2, we have xdC) = (^C) (1 + o(l)) as ^ + 00 . Hence we conclude that, as 

C —^ + 00 , 


V'2(a:c(C)) 


4 / 4 \ 

— log -C + 0 ( 1 ) 

a \oi } 

ilogC + 0(l). 
a 


Now we have the tools to estimate (3.13) and (3.14) for large values of C,. Dividing 
(3.13) by (C, taking liminf, and noting that when taking this limit all the (constant) 

C->+oo 

terms containing "Tt, after being divided by C, converge to zero, we obtain 


lim inf 


i(n — 1 ) 


> 


C^+00 V (n + 2 )CS?+" 2C,xl 

liminf i (C- T + f^(l-e)(logC- logT) - V'„(a 

C—>+oo L \ 


> 


(3.15) = 1 + lim 

C—>-+oo 

= 1 . 


Vn+2 D 


-e)MC 


-x{T+^{l-e)\ogT) 




Analogously, dividing (3.14) by C and taking lim sup we conclude that 

C->+oo 


lim sup I - _ 

>--1-00 \ (n + 2)^x, 




n(n — 1 ) 

~wr 


< 


(3.16) < 1 + hm^ (f^(l + e)i2^ - 1 (t + f^(l + e) logr) 




= 1 . 


Thus, (3.15)-(3.16) imply the existence of a function such that F’ri(C) = ^*( 1 ) 
for large C, and 


(3.17) 


(n + 2): 


Xn+2 


n(n - 1)\ 1 logC ^ 

- ^~2 - j ^ ^ - ~^niC)y as C —t + 00 . 


Let /3 := l/(n + 2), A := n{n — l)/2 and write (3.17) as 


(3.18) 


aP + Ax'^ = ( 1 + ^^^^A„((() ) , as ( + 00 . 
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Since —>■ 0 as C —>■ +oo we can use the geometric series expansion to write 

+ ' = + asC^+oo, 

and thus, from (3.18), we conclude that, as C —t +oo, 

(3.19) (l + ^x-) (l - ifiK(C) + h.o.t.) , 

where “h.o.t.” denotes terms with order higher than the orders of those explicitly 
written down. Taking into account (3.9) to estimate x" in the right-hand side of 

(3.19) , after some manipulations and using the binomial expansion we obtain, as 
C —t -boo, 

(3.20) ^ (t) "" “ + h.o.t. 


From standard center manifold theory [3, chap. 2], the long-time behaviour of 
(x(C),x(C)) is determined by the behaviour on the center manifold modulo exponen¬ 
tially decaying terms O (e“'''^) , where A € (0, a), in particular we can write 

^(C) = Xc(C) + O (e“^^) , 


and using (3.20) and remembering the definitions of /3 and A, we conclude that 


(3.21) 


x(C) = 


( “ 1 

1 , n{n - 1) , 

( “ ^ 

\(n -b 2)() 

' + 2a ' 

\in + 2)(J 


rt + l 
n+2 


+ h.o.t.. 


In order to obtain the corresponding estimates in the original time variable t we 
first need to relate the asymptotics of both time scales. 

Lemma 5. With the same notations and definitions as before, we have 


(3.22) 

and 


, _ (aPf .i-B 
^ ~ 1-/3 S 


as t ^ -boo, 


(3.23) C 



t^-f> -bo 



as C -boo. 


Proof. From (3.21) it follows that, for all e > 0, there exists T = T{e) such that, 
for all ( > T, 


(3.24) 




C J J J 

By the definition of f in (3.3) we can write 


C [1 — £, 1 -b e]. 


' L FA*' 
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and so dC,/dt = l/x{t). Hence dt/dC, = x{C,), and upon integration 

(3.25) t{C)-t{Co)=[ x{s)ds, 

JCo 

where Co = C(^o). Using the upper bound in (3.24) we have that, for Co > T, 

i(C)-i(Co) < (l + e)(a^)^”+^^'^- + (a/3)^ f^-f^ds 

^ J Co J Co 

= (1 + £)(a/3)-^H(C^ - Co^) + ^(C^-^ - Co-""), 

from which we conclude that 

^(C) < ^1% + (1 + £) + ^( 1 ) c oo. 


By considering the lower bound in (3.24) we would get the reversed inequality, with 
1 + e replaced by 1 — e, and since e is arbitrary, combining these two inequalities we 
conclude (3.22). 

In order to obtain C as a function of t we first consider the first term in the 
expansion (3.22), written in the form t(C) = (a/3)^(l — + o(l)), which 

implies C^“^ = — /3) t (1 + o(l)), and we substitute this into the lower order 

term C^ in (3.22), obtaining, as C —>■ oo, 


m 


.1-0 

{al3Y .1-0 
i-r 
{aPY .1-0 
i-r 

{aPY .1-0 


+ A{al3)-^Y+o[Y) 

+ A{aP)-P {Q^-P)^+o{Y) 

+ A{al3)~^ [{aj5)~^{l - f5)t{l + o(l)))’^ + o 


Write the last expression as 


= {I - P){aP)-H 


1 - A 


n + 1 


1 

n + 1 


n 

t "+1 + o 



Raising both sides to the power and using Newton’s binomial series in the right- 
hand side we obtain (3.23).D 

To finish the proof of Theorem 4 we use the expression (3.23) for ((t) in (3.21). 
After a few rearrangements similar to those above we obtain (3.1). □ 

4. Monomer long time behaviour in the modified time scale. As pointed 
out in Section 2, the proof of Theorem 2 relies in the exploration of (2.5), for which one 
needs detailed information on the long-time behaviour of ci(r). This will be obtained 
in this section, whose main result is the following: 

Theorem 6. With t and'ciir) as before, the following hold true: 


(4.1) (: 


Wn-l)/n 


a / 


(ci(r)) 


n —1 


l + (n-l)(l 


1) 

nJ J. 



as r ^ (X). 
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To prove this theorem we need to express t in terms of r in (3.1), and we shall 
do this in the next few lemmas. The approach is analogous to the one used in the 
previous section in order to relate t and C- 

Lemma 7. For T(t) defined by (2.2), and being ci(t) the first component of any 
solution (ci, Cn, Cn+i, ...) of (1.1), the following holds 

(4.2) = +2^logt + oi\ogt), ast^^. 


Proof. Let us consider the time scale (2.2) and write it as 

Cl(s)ds = r(T) + J ci{s)ds. 

Remembering that x{t) := ci{t), we deduce from (3.1) that, for all e > 0, there exists 
T = T{e) such that, for all t > T, 

*■*) Cl - 1 + ^1. 



Let us consider the upper bound case in (4.4) (the lower bound case is analogous). 
Substituting it in (4.3) one gets 


'(t) - t(T) = J ci{s)ds 


< 


a 


n + 1 


/■*_n(n—1) /■* 1 , 

/ s "+i(is + (l + e)-;— / -ds 

Jt n 1 Jrp s 


- T^) + (1 + log A, 


n+1 
n \ n+1 J 


from which follows 


r{t) < (+p) t+T + (1 + e)2^ \ogt + 0(1), as t ^ oo. 


By considering the lower bound in (4.4) we obtain the reverse inequality, with 1 + e 
substituted by 1 — e, and the result (4.2) follows by the arbitrariness of e. □ 

Lemma 8. With the assumptions of Lemma 7 we have, as r —>■ oo, 


t{r) = (+ t )"^ 

(4.5) 



(n - 1 )(+t)'^ logr). 


Proof Let B := (n^) ^ ^ write (4.2) as 

(4.6) r(t) = BL+i + Il(logt)(l + o(l)), ast,r—>-oo. 

Thus r(t) = Bt^^ (1 + o(l)), and so t = 5“+^(1 + o(l)). Substituting this 
expression in the logarithm term of (4.6) we get, as t, r —> oo, 

r(t) = Bt^^ + Dlog (^B~^ T~^ (1 + o(l))^ (1 + o(l)) 

= Bt^ + Dlog (^B~^T^~ ^ + o(logr). 
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and hence 


t "+i = i? (r — ^^Ulogr + o(logr)) ^ 

= Bt-i (l- (ll^i^r-ilogr + o(i2p)))”' 

(4.7) = Br-i (l + ^Dt-^ log r + o(i2^)) . 

Raising (4.7) to the power —and using Newton’s binomial series in the right-hand 
side we arrive at 


(4.8) t = B ^ r " — DB " t" logr-|- o(t" logr), 

and substituting B and D back into (4.8) results in (4.5). □ 

We are now ready to prove Theorem 6. 

Proof, (of Theorem 6). From (2.3) and (3.1) we have, as t,T —>• -|-oo, 



(=fi; 



i 

Vn-l-l 

n — 1 

\ n+1 n—1 

1 t "+i X 

X 

(> + 

n(n-l) /n-|-l\73 
n+1 \ a / 



n-l 


gF) 

n + 1 n-l n-l 

T " t 

X 


n(n—l)^ /n+l\n 

n+1 \ a / 


n—1 


n—1 


where we used Newton’s binomial expansion in the last equality. We can now apply 
Lemma 8 to write the right-hand side in terms of r, and after a few computations 
analogous to those described above, (4.1) arises. □ 

5. Rate of convergence to the similarity profile. In this section we prove 
the paper’s main result. Theorem 2. 

In [7] it was proved that the limit of (1.2) is equal to zero, thus generalizing a 
similar result first proved in [5] in the special case n = 2. For the proof of Theorem 2 
we shall apply the approach used in [5, 7], consisting in exploring the representation 
formula (2.5). In those papers the information that was needed was just the first 
term in Stirling’s expansion of the gamma function (to estimate the factorials) and 
the first term in the long time behaviour of the solution ci which, in [5] was obtained 
by center manifold methods, and in [7] by asymptotic analysis methods. Here we shall 
essentially make use of the results of Sections 3 and 4, in particular (4.1) proved in 
Theorem 6, in order to get the needed higher order information. 

Before starting the proof of Theorem 2 we briefly point out some of its conse¬ 
quences. 

Corollary 9. From Theorem 2 we immediately conclude the following: 

i. The convergence is not uniform in rj G (0, oo) \ {1}. 

ii. For 7] G 'E., a compact set not containing 1 and such that S fl (0,1) ^ 0 and 

S n (l,oo) 7 ^ 0, j/ 1 < /X < 2 — i, the solution converges to the similarity 
profile $1(77) at a power rate t~ ^. 
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iii. Under the same eonditions, ^ > 2 — ^ the solution converges to the similarity 
profile $i(? 7 ) at a rate 

Remark 3. Observe that Theorem 2 implies that for initial conditions (cj(0)) 
decaying fast enough with j, namely, for > 2 — the information about the initial 
data is lost in the similarity limit, in the sense that neither the bounds pi and p 2 , nor 
the decay rate p, and obviously no other details of the initial condition, are reflected 
in the limit or in the rate at which the limit is approached. 

However, if the initial data is slowly decaying, namely if l<p<2—then, 
although information about the initial condition is also lost in the similarity limit, 
the rate at which this limit is approached still retains information about the initial 
condition. 

Furthermore, if one can measure the decay at a single value of the similarity 
variable rj > 1 (and not only the overall rate in a set S of values of p) then the rate 
of decay of the initial condition for large clusters is always reflected in the rate of 
convergence to the similarity profile ^i{p). 

We think this is an interesting behaviour that, as far as we know, seems to be the 
first time it is observed in coagulation type systems. 

Proof, (of Theorem 2). As in [5, 7], we start by considering the case of monomeric 
initial conditions, in which case the sum in the right hand side of (2.5) is zero for all 
time T. 

Consider the function pi : [n,oo) x [0,oo) defined by 

/ nT_\U~U/^ pT 

(5.1) -.= (dflr-s)r-^s^-e-^ds. 

r(a;-n + l) Jq 

Comparing this expression with the integral in the right hand side of (2.5) we see 

that, when x = j £ Nn[n,oo), ^Pl{j,T) = (^)*'” ^^^"cj(T). So, we are going to study 
the rate of convergence of ipi{pT,T) to $i(??) as r —oo, for fixed p G (0,oo) \ {1}. 

Changing the integration variable in (5.1) s >—>■ y = s/t, and using (4.1) to write, 
for 0 < 2 / < 1, 


where /„(t (1 - y)) = (n - 1 ) (l - i) + o as t -)> oo, the expression 

(5.1) can be written as 


(5.2) 


>Pi{pt,t) 


.^l-n+7]T rl gT(t7logy-y) 

r(nr - n +1) X n+ /■("(! 


Now, in order to proceed, we need the following estimates, obtained from Stirling’s 
asymptotic formula, and valid as r —> oo, 

p, _ , IN _ r(?7T) _ 

” {riT - 1){VT - 2) ■ ■ ■ {pr - {n - 1)) 

" -V-' 

(n- 1 ) terms 

e-vr(^pr)vr-^^(l ^ + 0(t"2)) 
(^p)n-l(l_lhl^^+0(p-2)) 
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(5.3) 






1 - 


n(n—1) 1 


2 r}T 

Using the formula for the sum of a geometric series we can write 


0 {t 


-•2) 


1 - 


n 


(n - 1) 1 


+ 0{t-^)] 1 + 


(5.4) 


= 1 - 


Z r]T 
n{n —1)1 
2 rjT 

1 n{n — 1)\ 1 
12 


12r7T 


+ 0(r-2) 


-1 


+ 0(r-^) 1- 


12?7r 




rjT 


0{r 


- 2 \ 


Substituting (5.3) into (5.2) and using (5.4), equation (5.2) becomes 


ipi{r]T,T) = 


1 




1 _ 


1 + 6n(n —1)1 


X / (l+/n(T(l-y))) 


12 ryr 

gr(?7logy-y+7)) 




+ 0(r-2) ) X 

:dy 


(5.5) = Ji{r],T) - J 2 {v,t) + J3{v,t) - J4{v,t), 

where Jk{r],T), fc = 1,..., 4, are defined by 


(5.6) 

(5.7) 

(5.8) 

(5.9) 


Jiiv,T) := 

■= 

■= 

J4{v,t) := 


r{rjlogy-y+r]) 


y/2Trr] Jq ?/"'(1 — ?/)("' 


:dy, 


yn-yr ^ / J _ J) 


y/2TTr] 
X 


T2 


1277 T 


1 gT{ylogy-y+v) 


+ 0(r-") X 


-dy, 


'0 y-(l - 2 /)(—1)/" 

yri-rjT ^ pi griylogy-y+ri) 

_i f 1 + 6nin - 1) , 


■rfy, 


I 277 T 


\/2mj 

X [ fn{T{l-y)) 

Jo 


+ 0(r-") X 


^r(ylogy-y+y) 
1(1 _ y)(ii-l)/i 


:dy- 


We now study the limit as r —>• 00 of each of these functions separately. 
Lemma 10. For each fixed 77 C (0, 00 ) \ {!}, the following holds: 

Mv^t) = (4>i(77) + 0{1)t~^) 1(o,i)(?7) +ai{ri,T), 


as T ^ 00 , and ai{r],T) — >■ 0 beyond all orders. 

Proof. The proof of this result follows [5, Section 5.1] and [7, Propositions 6 
and 7], and is exactly the same as those proofs with the function if in those papers 
identically equal to 1. We need only to take the additional care of keeping track of 
the higher order terms, in order to be sure they are indeed negligible relative to the 
dominant terms arising from the other Jk in (5.5). 

If 77 > 1 and r > the function y 1 — >■ (77T — n) log y — yr is strictly increasing in 

y £ (0,1). Thus, for these values of r, and 
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we can estimate (5.6) as follows: 




pT(t7logy-y+i?) 


rT 2 




:dy 


1 


< „-TivioSV-V+i) 

- 

^ .jyn.^^g-'r(»7log?7-?7+l) 


f 


(1 - y)("-i)/’ 


:dy 


V27r?7 


If r? C (0,1) the approach is slightly more involved, but is also essentially the one 
used in [5, 7]. Fixing an e C (0,min {be \l-b}) , we decompose the integral in 
(5.6) as follows, 


(5.10) 




j^n-rjr ^ / re rl-e rl \ (? 7 log y-y+T?) 

72^^ Vo ^ Je 2/"(l - 


For the contributions arising from the first and third integrals in (5.10), the 
computations in [5, p. 385] show that both can be bounded above by a function 
which exponentially decreases to zero as r —> oo, namely: 


(5.11) 

(5.12) 


\/2Trrj 


+ 


'0 


gT(i?logy-y+i)) 

1-eJ y”(l - y)(”-i)/’ 


rdy = 


(1 _ ^lV)g"(l-log»)) ] 
n pl/^ 


Vy2^ 

/ nr\' 


+ 


\ yJ 2 ' KT \ (1 - e)*^ 


T-5g-TS3(l-e) 


for all sufficiently large r, where* goW) '■= iv^ogr] — rj) — (77 logy — y). 

For the remaining integral in (5.10) we again follow [5], as modified by [7]. With¬ 
out any change to what was done in those papers, an application of Laplace’s method 
[9, Chap. 3] results in the following, valid when r —>• 00 , 


(5.13) 


^n-r)T ^ rl-e ^r(r]logy-y+7]) ^ 

^/2'Kg 7g y"'(l — y)("“i)/" ^ (1 — 7y)("“i)/" ^ ^ 


This completes the proof of the lemma. □ 

Lemma 11. For each fixed rj C (0,oo) \ {1}, the following holds: 

Mv,t) = C>( 1 )t"^ 1 (o_i)( 77 ) -f 0:2(77,1-), 


as T ^ 00 , and 02 ( 17 , r) —>■ 0 beyond all orders. 

Proof. The proof is exactly equal to that of Lemma 10, with the exception that 
now all the terms will be multiplied by Hence the result imme¬ 

diately follows. □ 

Lemma 12. For each fixed 77 C (0,oo) \ {1}, the following holds: 


MV:r) = (tz- 1 ) (1 - i) 


log((l - 77 )r) 


(1 — 77)("“^)/"' (1 — 77)t 


(1 -Lo(l))l(o,17(77) -F 03(77, r), 


This notation was employed in [5, 7] and is maintained here for convenience of the reader. 



















16 


F. P. DA COSTA, J. T. PINTO AND R. SASPORTES 


as T ^ oo, and a^{r],T) —>■ 0 beyond all orders. 

Proof. The proof of this lemma follows the steps taken in the proof of Lemma 10 
with the function inside the integral now containing the additional multiplicative 
factor 


/n((l - y)T) = (n - 1) (1 - i) + o , as T oo. 

As before, the cases rj > 1 and rj G (0,1) need to be treated separately, and the 
last one requires the splitting of the integral into three terms, two of them dealing 
with regions near the boundary. Since 0 as t —oo, the asymptotically 

exponentially small bounds when t —>■ oo obtained in the proof of Lemma 10 for the 
cases ?7 > 1, and for the “boundary integrals” and in the case ry G (0,1), still 
hold in the present case. 

Thus, we are left with estimating the integral which can again be done 

by Laplace’s method exactly as before, giving the same result as (5.13) with the 
additional multiplicative factor This completes the proof. □ 

Lemma 13. For each fixed rj G (0,oo) \ {!}, the following holds: 

Mv,t) = ^l(o,i)(»7) + a4(»7,T), 

as T ^ oo, and 04 ( 77 , r) —>■ 0 beyond all orders. 

Proof. The proof of this lemma is related to that of Lemma 12, as the proof of 
Lemma 11 is to that of Lemma 10: the expressions of the corresponding Jk differ by 
the factor which entails the result. □ 

Putting together the results in Lemmas 10-13 we conclude the proof of Theorem 2 
in the case of monomeric initial data, that is: when Cj(0) = 0 for all j > n, and the 
only initial component potentially nonzero is ci(0). 

Let us now consider the case when the initial data is nonzero. More specifically 
we shall consider initial conditions satisfying the following assumption: 

(HI) 3pi,p2 > 0, fi>l : Vj > n, j^Cj{0) G [pi,P 2 ]- 

So, let us now consider the sum in the right hand side of (2.5). We are interested 
in studying the behaviour of 


V a 


g-l J 


E 

k—n 


rJ-k 


ij-ky. 


Cfc(O), 


when j, T —> + 00 , with rj = jfr 1 fixed. As in [5, 7], we shall write v := rj~^ 
and use t = ju in the above expression. We are then left with studying the limit, as 
j —>■ 00 of 


(5.14) 



k—n 


(j - ky 


Cfe(O). 


Consider the cases tz > 1 and v G (0,1) separately. 

When 7z > 1 the limit is easy to handle, giving an asymptotically exponentially 
small contribution. In fact, the situation is entirely analogous to the one in [5, page 
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387]: using the upper bound in (HI) and changing the summation variable fc i—>• £ := 
j — k in (5.14) we get 


,5.15, 



k—n 


(j - k)\ 


Cfe(O) < 


P2_ 



j n e ^ ] 




£! 


and from here the computations in [5] apply almost verbatim since (5.15) with n = 2 
reduces to the expression in [5]. The hnal result is the bound 



k—n 


U - k)'. 


Cfc(O) < 




e-^('^-i-^°g‘^)(l + o(l)), 


(5.16) 

as j —>• oo. Since > 1, we have v — 1 — iogiy > 0, and thus the upper bound is 
asymptotically exponentially small in the large j limit. 

When 1 / G (0,1) we again follow [5, page 388]: again changing the summation 
variable k i := j — k and using the upper bound in (HI), we can choose a 
fixed P G min{i/e, 1}) and split the sum into a “small” £ and a “large” £ 

contribution. 


E 


£=0 


£!(j - £Y 


^ , y- {J^Y 

^ £\(j-£Y^ ^ £\{j-£Y' 

0<£</3j ' pj<l<j-n ’ 


The analysis of each of these terms proceeds exactly as in [5, page 388], and concluding 
that the contribution from the “large” I is asymptotically exponentially small when 
j —>■ oo, and the one due to the “small” I sum being bounded above by 


(5.17) 


P2 


((1 — P)ny \a 




where we have used the fact that jv = t. 

To complete the proof we need to obtain a lower bound for the limit of (5.14) 
when j —>■ oo. Using the lower bound in (HI) and the same change of variable fc i—>• £ 
as above, it is easy to obtain the bound 


,5.15, (!^)- 


Li •? 


{jB 


/ • \n _ Ic n — 1 3 — 71 

k=n 


Now write 


3 \i ^ 


e=o 


e=o 


£! 


E 


i-j—n-Yl 


BY 

£! 


Since i' G (0,1), the first term in the right hand side converges to 1 as j —oo (cf., 
e.g., [1, formula 6.5.34]). For the second term, the sum of which has n terms, a finite 
and fixed number independent of j, we have 


3 

^ E ■ 

n+l 


BpY 

£! 






(i-i)! 


(i-2)! + + O-n+l)! ) 
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+ i + ^ + + (i + o(r-)) 

*^(ilSS-r(i + o(r‘)). 


Using Stirling’s expansion we have e ^ 


-iv {jvy _ _ 1 


.eJ(i-‘'+i°g'^)(l + e>(j-i)). From 


]\ - 

•3 (juY 

V S (0,1) it follows that 1 — v -\- \ogv < 0, and thus e~^'' ^ —-— converges to 




£! 


zero exponentially fast as j —> oo. This implies that (5.18) can be written as 


(5.19) (( 






k—n 


\ a J 


Hence, a lower bound for the limit of (5.14) as j > oo is 


(5.20) pi0) "r?-'^r^-^(l + o(l)), 

where we have used the fact that jV = r. From (5.18) and (5.19) the term (1.4) 
follows. 

This concludes the proof of Theorem 2. □ 
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